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Abstract 

Recent exact predictions for the massive scaling limit of the two 
dimensional XY-model are based on the equivalence with the sine- 
Gordon theory and include detailed results on the finite size behavior. 
The so-called step-scaling function of the mass gap is simulated with 
very high precision and found consistent with analytic results in the 
continuum limit. To come to this conclusion, an also predicted form 
of a logarithmic decay of lattice artifacts was essential to use for the 
extrapolation. 
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1 Introduction 



For two-dimensional quantum field theories there are exact solutions avail- 
able in quite a number of cases. The meaning of the notion 'exact' is how- 
ever not quite universal. In many cases it does not imply a straight-forward 
evaluation, and often not rigorously provable conjectures are adopted at in- 
termediate steps. To corroborate such solutions it is very attractive to com- 
pare — whenever possible — high precision numerical simulations with exact 
predictions. An interesting such opportunity arose in recent years from an 
exact solution of the famous XY-model |2j in its massive continuum limit 
taken from the vortex-phase. This solution involves a conjectured chain of 
equivalences that links it to the exactly solved sine-Gordon field theory 
The predictions are particularly suitable for a numerical check since they 
include information about the XY-model in a finite volume and even about 
the asymptotic lattice spacing dependence for the standard discretization of 
the model. On the simulation side the XY-model profits from the practi- 
cal absence of critical slowing down for cluster algorithms [3], allowing for 
large correlation lengths to be simulated, and from the availability of reduced 
variance estimators for correlations [5 . 

Interesting light is shed on a problem inherent in all numerical lattice 
field theory computations, including QCD and asymptotically free two-di- 
mensional models: the need of a continuum extrapolation. To carry it out, 
a theory based analytic form for the asymptotic dependence on the lattice 
spacing a is required. Usually a power behavior a and/or a 2 from Symanzik 
theory jH] , based on all-order perturbation theory, is used. In [3 |H] high pre- 
cision simulations of the 0(3)-model revealed however that the expected a 2 
behavior does not prevail down to very small lattice spacings, and at present 
there seems to be no theoretical understanding for this. The XY-model is 
not asymptotically free in the usual sense and the exact prediction for its 
a-dependence differs from the Symanzik behavior, which is very interesting 
to verify. A good understanding of the cut-off dependence is clearly impor- 
tant also for QCD simulations, in particular now that they become more and 
more precise and hence quantitatively useful for phenomenology. 



2 



2 Theoretical predictions 



2.1 The LWW-coupling 

The two dimensional nonlinear 0(2) symmetric cx-model can be defined by 
the classical Lagrangian 

£=^ M S ,a fyS ,a , S a S a = l, a = 1,2. (1) 
2g l 

When space-time is discretized the XY-model is recovered, its standard ac- 
tion being 

S = -/3^2co8(9 k -6 l ). (2) 

{k,l) 

The summation is performed over all nearest neighbor pairs. The model 
describes a set of spins with unit length arranged on the sites of a square 
lattice that are subject to a ferromagnetic interaction. The angle 9 n denotes 
the alignment of the spin at site n relative to an arbitrary axis and f3 = 1/g 2 
is the inverse temperature or bare coupling. 

Let the spatial extent of the lattice be L. This introduces a natural 
external scale and renormalized couplings that run with L can be defined. 
In |H] the Liischer-Weisz- Wolff-coupling (LWW) was introduced 

g\L) = 2M{L)L. (3) 

Here M(L) is the finite volume mass gap of the theory defined by the transfer 
matrix in the sector of zero spatial momentum. The running of the coupling 
is described by the /9-function 1 

( 4) 

For finite scale changes the universal step-scaling function a has been defined, 
and it describes how the LWW coupling changes if the scale is dilated by a 
factor s 

a(s,f(L))=f(sL). (5) 

The lattice version (a is the lattice spacing) of the step scaling function 
T*(s,g 2 ,a/L) depends in addition on the resolution of the lattice. 

lr The symbol (3 is traditionally used for both the inverse temperature and the j3- function 
which should not be confused with each other. 
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2.2 DdV equation 

Starting from the light-cone lattice regularization ^01 E] of the sine-Gordon 
model and the corresponding Bethe Ansatz solution, a complete descrip- 
tion of the exact spectrum of the model in finite volume is provided by the 
Destri-deVega (DdV) non-linear integral equations. They were originally 
suggested O H3| for the ground state of the sine-Gordon or the closely 
related massive Thirring model, but can be systematically generalized to 
describe excited states as well; the (common) even charge sector of both 
models [T^l CS1 HE] and the (different) odd charge sector in both sine-Gordon 
and massive Thirring model [T7] . 

The ground state energy is given by 

+00 

Eq{L) = -— Im J dx sinh(x + vq) In [l + e iZo{x+t ^] . (6) 
—00 

The mass M without an argument denotes here the infinite volume mass gap 
M = M(L = 00). The parameter rj is needed to keep the branch of the 
logarithm well-defined. The result is independent of rj as long as < r] < n. 
The so-called counting function Z (9) solves the nonlinear integral equation 

Z Q (9) = ML sinh 9 



— 1 



/ <MG (,-,- i ,)- G ( J+ , +j ,) 1M i + ^] 

—00 

with a kernel given by 

, , 1 7 ho sinh^^ 

G{9) = — dk e lke ^ ^ , 

2tt J 2 cosh 4 sinh ^ 



(7) 



where p = P^ G /(8rr — j3^ G ) parameterizes the sine-Gordon coupling /3 s g. 
Similarly the energy of the first excited state is given by 

+00 

E 1 (L)=M-—lm J dx smh(x + ir]) In [l - e iZl(x+iv) ] , 



(9) 



with Z\(9) being a solution of 

Z x {9) = ML sinh 9 + x(0) 
+00 

-i J dx [G(9 -x-ir))-G{9 + x + ir])} In [l - e iZl( - x+iv) ] . ^ 
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The source term is of the form 



X(6) = 2tt dzG(z) 



(11) 



and r) is arbitrary within the same range as before. 

The finite volume mass gap of the theory is M(L) = E\(L) — E (L). The 
XY-model is believed to lie in the same universality class as the sine-Gordon 
model with the sine-Gordon coupling /3 s g = VSn, or p — > oo [3J. In this case 
the kernel simplifies to 



G{9) 



1 

2^ 



AkO 



dk 



e 7T[fe| _|_ I ■ 



'12) 



The Fourier integral can be carried out 



G{9) 



Air 2 



T .n-i9. T .ir + i9. T ,„ i9 . _ _ . 



(13) 



in terms of the digamma function 2 ^(z). Also the source term can be calcu- 
lated exactly in the limit p — » oo, 



x(x + z?y) = — i In 
+ i In 



1 _ -r] 

2 2tt 



r i 



IX — T] 

2vr 



+ z In 
— i In 



1 — 77 
1 2 2tt 



r 1 + 



IX — T) 

2vr 



(14) 



2.3 Asymptotic behavior 

For large volumes (L — > 00) and also for small volumes (L — > 0) asymptotic 
formulae exist which can be used to calculate the finite volume mass gap in 
this model independently of the DdV equation. 

The ftsQ — > 8ir sine-Gordon model is equivalent to the SU(2) symmetric 
chiral Gross-Neveu model and is thus asymptotically free. In principle it is 
straightforward to do perturbation theory for the XY-model using the for- 
malism of Amit et al. The final result, as usual in an asymptotically free 
model, is a power series in a running coupling. Technically however, the cal- 
culations are more involved than usual, since the perturbative contributions 
cannot be written in terms of usual Feynman integrals. 

2 The digamma function Vf'(z) is the logarithmic derivative of the gamma function, given 

by *(z) = r(z)/r(z). 
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The 2-loop perturbative result for the LWW coupling is 

^ = ^(l + 7 + W) + 0(7 3 ), (15) 
with the running coupling 7 solving 

— - ~ln(2 7 ) = In ( — - ) . (16) 



2 7 2 v " \ML 

Along the lines described in JH] we obtain 

u 2 = 0.20(1), (17) 

where due to the complexity of the computation the 2-loop coefficient has to 
be obtained numerically in the end. 

On the other hand, Liischer has found a universal formula which describes 
the leading large volume corrections to the mass of the lightest particle on 
a spatial torus in terms of scattering amplitudes of the theory The 
formula exists in any dimension. It is particularly useful in two dimensional 
integrable models where the scattering amplitudes are known exactly. In our 
case it gives 



tf=2ML + — IdB cosh£e- MLcoshe 



+ <D(e 



71 

-2ML\ 



1 + e 



ix(e+ ! f) 



3 Numerical results 



3.1 Numerical solution of the DdV equation 

To calculate the finite volume mass gap M(L) it is necessary to evaluate the 
nonlinear integral equations and ffTTTf) numerically. The idea is to solve 
them iteratively [20] • In a first step the integral term on the right hand side 
is neglected, in every consecutive step the previous approximation is inserted 
into the integral. The iteration is aborted when the desired precision is 
reached. Further details of this procedure are described in [21]. Figure ^ 
shows some values of the LWW-coupling as a function of the the size in units 
of the infinite volume mass gap obtained from the numerical solutions of the 
DdV equations and compares them to the asymptotic formulas (JT3J) and (|TS|) . 

To compute the step-scaling function a(2,u) first a value ML has to be 
found such that after solving the DdV equation for it 2M(L)L = u holds 
accurately enough. Then a is given by solving again with the doubled value 
2ML. 
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Figure 1: LWW coupling obtained from solutions of the DdV equation. 



3.2 Monte Carlo study of the model 

The XY-model can be simulated very effectively with cluster algorithms [H 
E21 1231 • The finite volume mass gap may be obtained from the decay of the 
correlation function of spatially averaged spin fields (see 9 J for details). For 
our simulations we use a single cluster algorithm and an improved estimator 
for the correlation function j^]. 

To calculate the step scaling function we perform calculations on lattices 
of different resolutions a/L. For each lattice a bare coupling (3 is determined 
that leads to the desired value u of the LWW-coupling. A simulation at the 
same bare coupling but with lattice size 2L/a yields a point of the lattice step- 
scaling function E(2, u, a/L). To obtain the continuum step-scaling function 
it is necessary to perform a continuum extrapolation. The unusual form of 
the lattice artifacts in the XY-model was predicted in [23 and in the case of 
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the step-scaling function reads 

£(2, g\ a/L) = a(2, f) + ° + + 0(ln" 4 , (19) 

where U is a non-universal constant (U = 1.3(1) for the standard action [2*K] ) 
and c can be calculated for each g 2 (see appendix EJ. The lattice resolution 
is expressed in terms of the infinite volume correlation length £ = = j^jr, 
which can be measured directly if it is not too large. When the bare coupling 
P approaches f3 c , the Kosterlitz Thouless (KT) theory |IJ E] predicts the 
correlation length to diverge according to 



£ = Aexp [C 



Pc~P 



Pc 



(20) 



The non-universal constants are known with good precision for the standard 
action [2l)] . 

A = 0.233 ±0.003 (21) 
C = 1.776 ±0.004 (22) 
p c = 1.1199 ±0.0001. (23) 

This result was derived in an approximative renormalization group study and 
is strictly valid only for P — > P c . At P ~ 1.05 we see a numerical matching 
with KT, so we use eq. f!20|) to determine the correlation lengths yet closer to 
criticality. We emphasize that this result is only used for the lattice artifacts 
and that errors in £ are suppressed by the logarithm in (fT9|l . 

Fitting the Monte Carlo estimates for the lattice step-scaling function to 
(JEJ) yields a point of the continuum step scaling function and the constant c. 
Both values can be compared to predicted values, which is done for several 
values of the LWW-coupling in table ^ The small discrepancy between 
Cth and cmc can presumably be attributed to subleading cutoff effects in the 
Monte Carlo data. Figure El illustrates the extrapolation for one of the points, 
the corresponding figures for the other points look qualitatively the same. 
Figure El which has an equally good x 2 - value, demonstrates how important 
the knowledge of the lattice artifacts is. The data needed to extract these 
results are collected in table El of appendix |B] 



4 Conclusions 

In this work we have found a remarkable consistency between analytic and 
numerical results. For several values of the LWW coupling we found per- 
fect agreement between numerical continuum-extrapolated values of the step 
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t 


^Ddv(2,^ 2 ) a MC (2,g 2 ) 


Cth CMC 


X 2 /dof 


3.0038 
1.7865 
1.6464 
1.6020 


4.3895 4.40 ±0.02 
1.8282 1.829 ±0.007 
1.6515 1.657 ±0.003 
1.6029 1.608 ±0.004 


2.6176 ±0.0002 2.4 ± 0.6 
5.30 ±0.01 4.8 ±0.5 

5.4 ±0.2 4.3 ±0.3 

5.5 ±1.5 4.4 ±0.5 


2.51/3 
0.73/3 
0.35/3 
0.90/3 



Table 1: At different values of the LWW-coupling theoretical predictions 
(DdV) for the step-scaling function are compared with numerical results 
(MC). Also the constant c as predicted by theory (th) is compared to the 
numerical value. 



scaling function and those following from solving the DdV equation. This 
agreement is only obtained by employing for the extrapolation a logarithmic 
dependence on the lattice spacing that follows from the connection between 
the sine-Gordon and the XY-model together with the KT-behavior of the 
correlation length close to the critical point. 

This overall consistency furnishes strong evidence that all the considera- 
tions and assumptions leading from the six- vertex model via Bethe Ansatz to 
a set of integral equations that can be used to calculate the energy spectrum 
of the sine-Gordon model are correct. Simultaneously it also corroborates 
the assertion that the massive continuum limit of the lattice XY-model can 
be described by the 0l G — > 8tt limit of the continuum sine-Gordon model. In 
Ref. the latter was investigated by the Form Factor Bootstrap construc- 
tion. The renormalized 4-point coupling g-^ and 2-point correlation functions 
were compared to their lattice counterparts. Similarly to our findings here, 
after (and only after) the logarithmic lattice artifacts were taken into ac- 
count, non-trivial agreement between lattice data and Form Factor results 
was found. 

Indirectly also the KT-scenario is confirmed once more. Calculations lead- 
ing to this scenario are more or less approximate and there has been a con- 
troversy whether they can be trusted or not for a long time j22j . This result 
demonstrates in an impressive way, how important it is to have theory-based 
information on the functional form of lattice artifacts in order to extract the 
right continuum values from numerical data. Although the accessible range 
of lattices is much bigger and the precision higher than it would be in four 
dimensional theories, an extrapolation to the continuum could lead to un- 
controlled systematic errors for precise data, if the right form of the lattice 
artifacts were not known. 

Acknowledgement. We would like to thank Burkhard Bunk for discus- 
sions. This investigation was supported in part by the Hungarian National 
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Figure 2: Comparison of MC-data with a numerical solution of the DdV 
equation at g 2 = 1.7865. The spatial extents L/a of the lattices were 
10, 20, 40, 80, 120 and 160. The smallest lattice was discarded in the 
fit. 
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A Lattice artifacts of the step-scaling func- 
tion 

The LWW-coupling g 2 = 2M(L)L on a lattice is affected by lattice artifacts. 
In contrast to the simulations we now consider families of lattice calculations 
at fixed physical volume in units of the infinite volume mass gap, that is 
fixed ML. Then each chosen resolution a/L in principle fixes a value for (3 
and hence also g 2 . In this sense we find g 2 = /(ML,(), where £ = (Ma) -1 
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• MC-data 

* solution of DdV equation 
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Figure 3: Comparison of MC-data with a numerical solution of the DdV 
equation at g 2 = 1.7865. The predicted continuum value could not have been 
confirmed by an extrapolation with a second degree polynomial (including a 
linear term), as used in this plot. 



is used. The leading ^-dependence is predicted by theory [24J to be 



f(ML,0 = f (ML) - ' + 0([ln£]- 4 ), f (ML) = f(ML, oo) 



fi{ML)n 2 

(24) 

where U is an action dependent constant (U ~ 1.3 for the standard action) 
and fi(ML) a coefficient. In solving the DdV equation (see figure HJ) we 
evaluated fo(ML). In [25 it was shown that the same function f\ also char- 
acterizes the approach of the the sine-Gordon model to its XY-model limit. 
Here one may also define an LWW coupling g 2 G and study its dependence on 
p for fixed ML to derive 

f sG = f(ML,p) = f (ML) + 2/1 ^f L) + 0(p-') (25) 
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for large p. This relation is entirely in the continuum and it hence becomes 
clear that fi is universal. 

Numerically we avoid reference to infinite volume quantities (except here 
for the artifacts) and determine the step-scaling function a(2,g 2 ) at a certain 
fixed g 2 = u. The relation to the present language is 

a(2,u) = f [2f 1 (u)], (26) 

where f^ 1 is the inverse function of fo. At finite £ a computable slightly 
different value (M + 5M)L is associated with the target u- value. Then the 
doubled value is mapped to £ instead of cr. Collecting now all these correc- 
tions to leading order we obtain 



£(2,u,a/L) - (t(2,u) 



n 2 



2 



V ti(ML) 7 2(lnf + C7) 



which defines c(ML). 

Further corrections to (|2*7j) are of order 0([ln£] -4 ). The required deriva- 
tives of fo in the continuum can be obtained numerically from solutions of 
the DdV equation. Also the coefficient f\ may be gained from a numerical 
solution of the sine-Gordon version of the DdV equations. Instead of using 
the XY-model specific integral kernel the version (jHJ) with finite p is 
taken. All data needed to calculate c is collected in table El 
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B Tables 



In the following two tables some of the data is collected that is needed to 
extract the results. An estimate of the derivative Jp is used to propagate 
into the total error of £ the statistical error of g 2 and to correct the value 
of £ for the small difference between u and g 2 . In most cases this correction 
amounts to a small shift compared with the statistical errors. The results in 
table El were obtained twice independently. 



L/a 


P 




g 2 m 


u 


E(2, u, a/L) 


20 


0.9649 


3.0022(45) 


4.5144(74) 


3.0038 


4.519(13) 


40 


1.0000 


3.0036(41) 


4.4824(75) 


3.0038 


4.483(12) 


80 


1.0241 


3 0146(39) 


4 5055(76) 


3 0038 


4 480(12) 


120 


1.0359 


2.9919(39) 


4.4232(64) 


3.0038 


4.451(11) 


160 


1.0423 


3.0069(38) 


4.4742(64) 


3.0038 


4.466(11) 


10 


1.0358 


1.7866(16) 


1.9379(19) 


1.7865 


1.9378(31) 


20 


1.0541 


1.7837(16) 


1.9161(19) 


1.7865 


1.9222(40) 


40 


1.0668 


1.7860 16 


1.9010(19) 


1.7865 


1.9019 35 


80 


1.0756 


1.7896(16) 


1.8973(19) 


1.7865 


1.8911(37) 


120 


1.0800 


1.7882(17) 


1.8898(19) 


1.7865 


1.8869(38) 


160 


1.0830 


1.7850(16) 


1.8778(13) 


1.7865 


1.8796(21) 


20 


1.0758 


1.64585(73) 


1.71173(79) 


1.6464 


1.7125(12) 


40 


1.0844 


1.64632(75) 


1.70116(83) 


1.6464 


1.7013(16) 


80 


1.09076 


1.64646(73) 


1.69317(81) 


1.6464 


1.6931(13) 


120 


1.0937 


1.64636(74) 


1.68885(79) 


1.6464 


1.6889(13) 


160 


1.09551 


1.64645(77) 


1.68774(80) 


1.6464 


1.6877(13) 


10 


1.0720 


1.5971(12) 


1.6721(66) 


1.602 


1.6801(69) 


20 


1.0825 


1.6033(12) 


1.6574(13) 


1.602 


1.6556(21) 


40 


1.0903 


1.6032(12) 


1.6463(14) 


1.602 


1.6449(20) 


80 


1.0960 


1.6020(11) 


1.6397(13) 


1.602 


1.6397(20) 


120 


1.0984 


1.6034(13) 


1.6365(14) 


1.602 


1.6349(20) 


160 


1.1000 


1.6028(13) 


1.6330(13) 


1.602 


1.6321(20) 



Table 2: Simulation parameters and results. In the last column all sources 
of errors have been combined. 
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ML 


V 




fi 


a ml 


2.164 x 1(T 12 


100 
300 


1.60461349 
1.60229497 








oo 


1.60200011 


13.1(2) 


5.616(4) x 10 8 


4.328 x 1(T 12 


100 
300 


1.60541070 
1 60315200 








oo 


1.60286501 


12.7(2) 


2.962(4) x 10 8 


9.11895 x 10" 6 


100 
300 


1.64750389 
1 64652300 








oo 


1.64640000 


5.52(2) 


7.5826(6) xlO 2 


1.82379 x 10" 5 


100 
300 


1.65254276 
1.65162280 








oo 


1.65150749 


5.18(2) 


4.3054(4) xlO 2 


0.014024 


100 

9nn 
zuu 

300 


1.78689383 

1 78fi^QQzL^ 

1 78654489 








oo 


1.78650122 


1.963(1) 


3.6430(4) 


0.028048 


100 
9nn 

300 


1.82853577 
1.82824321 








oo 


1.82820661 


1.646(1) 


2.5230(2) 


0.9652 


100 
9nn 

300 


3.00389085 
3.00384633 








oo 


3.00384076 


0.25041(1) 


1.27592(2) 


1.9304 


100 
200 
300 


4.38947652 
4.38946271 
4.38946015 








oo 


4.38945811 


0.092061(3) 


1.58593(6) 



Table 3: Solutions of the DdV equation needed to calculate c t \ x . Numerical 
values of g 2 G are uncertain by one in the last digit. 
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